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Abstract. The autocorrelation function and similarly the Fourier-power spectrum of a rotation measure (RM) map of an 
extended background radio source can be used to measure components of the magnetic autocorrelation and power- spectrum 
tensor within a foreground Faraday screen. It is possible to reconstruct the full non-helical part of this tensor in the case 
of an isotropic magnetic field distribution statistics. The helical part is only accessible with additional information; e.g. the 
knowledge that the fields are force-free. The magnetic field strength, energy spectrum and autocorrelation length A B can be 
obtained from the non-helical part alone. We demonstrate that A B can differ substantially from Arm, the observationally easily 
accessible autocorrelation length of an RM map. In typical astrophysical situation A RM > A B . Any RM study, which does not 
take this distinction into account, likely underestimates the magnetic field strength. For power-law magnetic power spectra, and 
for patchy magnetic field configurations the central RM autocorrelation function is shown to have characteristic asymptotic 
shapes. Ways to constrain the volume filling factor of a patchy field distribution are discussed. We discuss strategies to analyse 
observational data, taking into account - with the help of a window function - the limited extent of the polarised radio source, 
the spatial distribution of the electron density and average magnetic energy density in the screen, and allowing for noise reducing 
data weighting. We briefly discuss the effects of possible observational artefacts, and strategies to avoid them. 

Key words. Magnetic Fields - Radiation mechanism: non-thermal - Galaxies: active - Intergalactic medium - Galaxies: 
cluster: general - Radio continuum: general 



1. Introduction 

1.1. Cosmic magnetic fields 

The interstellar and intergalactic plasma is magnetised. The 
origin of the magnetic fields is partly a mystery, yet it allows 
fascinating insights into dynamical processes in the Universe. 
Magnetic fields are an important constituent of cosmic plasma 
in so far as they couple the often collisionless charged particles 
by the Lorentz-force. They are able to inhibit transport pro- 
cesses like heat conduction, spatial mixing of gas, and propa- 
gation of cosmic rays. They are essential for the acceleration 
of cosmic rays. They mediate forces through their tension and 
pressure, giving the plasma additional macroscopic degrees of 
freedom in terms of Alfvenic and magnetosonic waves. They 
allow distant cosmic ray electron populations to be observed 
by magneto-curvature (synchrotron) radiation. 

Observational studies of spiral galaxies have revealed 
highly organised magnetic field configurations, often in align- 
ment with the optical spiral arms. These magnetic fields are 
believed to be generated and shaped by the dynamo action of 
the differentially rotating galaxy disks from some initial weak 
seed fields. 

The seed fields could have many origins, ranging from out- 
flows from stars and active galactic nuclei over battery effects 
in shock waves, in ionisation fronts, and in neutral gas-plasma 
interactions, up to being primordially generated in high en- 



ergy processes like phase transitions or inflation during the very 
early Universe. 

In order to learn more about the magnetic field origin, 
less processed plasma has to be studied. There is the possi- 
bility that the magnetic fields outside galaxies, in galaxy clus- 
ters and - if existing - even in the wider intergalactic space 
carry more information on the field's origin. In clusters, mag- 
netic fields with a much lower degree of ordering, compared 
to the organised fields in spiral galaxies, have been detected. 
However, they may be highly processed by turbulent gas flows 
driven by galaxy cluster mergers, which may mask their origin. 
Regardless, cluster magnetic fields are an interesting labora- 
tory to study magneto-hydrodynamical (MHD) turbulence, and 
are of great importance to understand thermal and non-thermal 
phenomena in the intra-cluster medium. 

Despite their obvious importance for many astrophysi- 
cal questions, and despite many observational efforts to mea- 
sure their properties, our knowledge of galactic and inter- 
galactic magnetic fields is still poor. For an overview on the 
present observatio nal and theo r etical knowledge the excellen t 
review articles by Rees ( 1987), Wielebinski & Krause ill 9931) 
iKronberd ill 994. iBeck et al] ill 9961) iKulsrudl J 1999b iBeckl 
bOQlh. Idrasso & RubinsternT^OOll) . ICarilli & TavloTj2002t) . 
and Widrow (2002) should be consulted. 
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1.2. Faraday rotation 

One way to probe magnetic fields is to use the Faraday rotation 
effect. Linearly polarised radio emission experiences a rotation 
of the polarisation plane when it transverses a plasma with a 
non-zero magnetic field component along its propagation di- 
rection. If the Faraday active medium is external to the source, 
a wavelength-square dependence of the polarisation angle mea- 
sured can be observed and used to obtain the RM, which is the 
proportionality constant of this dependence. Such situations are 
realised in nature in cases where a polarised radio galaxy is lo- 
cated behind the magnetised medium of a galaxy, or behind or 
embedded in a galaxy cluster. 

The focus of this work is on the analysis of RM maps of 
Faraday screens, in which the fields are statistically isotrop- 
ically distributed. This should be approximately fulfilled in 
galaxy clusters, but not in the highly organised spiral galax- 
ies. However, our analysis should also give some insight into 
the statistics of RM maps of galaxies, since many of the results 
do not strictly require perfect isotropy. 

Magnetic fields in galaxy clusters are know n to exist due t o 
detection of cluster wide synchrotron emission (Willson 1970), 
and detection of their Faraday rotation effect. Although the as- 
sociation of the RM with the intra-cluster medium is not unam- 
biguous, since it could also be prod uced in a magnetised plasma 
skin of the observed radio galaxy dBicknell et aljfl990l) . there 
are arguments in favour of such an interpretation: (i) The asym- 
metric depolarisation of double radio lobes embedded in galaxy 
clusters can be understood as resulting from a difference in the 
Faraday depth of the two lobes (IT ,ainJ 1 988l |p arrington lr"aTl 
1988). (ii) A recent RM study bv lClarke et alJ d200lh of point 
sources located mostly behind (but 40% inside) galaxy clusters 
show a larger dispersion in RM values than a reference sample 
without a galaxy cluster intersecting the line-of-sight. (iii) The 
cluster- wide radio ha los observed in some clusters of galax- 
ies fe.g. lFerettilll999h show synchrotron emission of relativis- 
tic electrons within magnetic fields. The cluster fields strength 
should be within an order of magnitude of their Faraday rota- 
tion estimates for the radio emitting electrons to have a reason- 
able energy density (compa red to the thermal one) (e.g. as can 
be read off Fig. 1 in EnBlin & Biermannll998l) . 

Typical RM values of galaxy clusters are of the order of a 
few 100 rad/m 2 , being consistent with field strengths of a few 
/iG which are well below equipartition with the thermal cluster 
gas. However, in cooling flow clusters extreme RM values o f 
a few 1000 rad/m 2 were detected (see Carilli & Tavlon l2002). 
indicating possibly substantial magnetic pressure support of the 
intra-cluster gas there. 

Although the magnetic fields of galaxy clusters are less 
ordered than these of spiral galaxies, the presence of co- 
herent structures is suggested by high resolu tion Far aday 
, which exhibit sometimes RM bands (e.g. iDreher et al 



maps, wmcn exhibit sometimes KM bands (e.g. Drener et al. 
1987l lTavlor fc Perlevll993UTavlor et all200lHFilek & Owen 



2002). Such bands may be caused by shear-amplification of 
originally small-scale magnetic fields, as seen i n numerical 
MHD simulations of galaxy cluster formation dPolag et alJ 
119991) . They are likely embedded within a magnetic power- 
spectrum which extends over several orders of magnitude in 



wavevector space. Similar to hydrodynamical turbulence, a 
broad energy injection range is followed by a power-law spec- 
trum at larger wavevectors. For attempts to measure the mag- 
netic £oror_srecttum from cluste r simulations and radio maps 
see lPolag et alJ J2002I) and lGovoni et alJ d2002l) respectively 1 . 

Another area of application of the theory developed here 
can be to measure the properties of an hypothetical large-scale 
magnetic field outside clusters of galaxies, which could be of 
primordial origin. A pioneering feasibility study in this direc- 
tion was done by KolatJ ll 19981) . who already outlined several 
of the ideas investigated in this work. He proposed to probe the 
cosmological magnetic fields by using catalogues of RM mea- 
surements of distant radio galaxies and to measure spatial RM 
correlations between them in order to measure the magnetic 
power spectrum, as we propose to do for extended Faraday 
rotat ion maps. A RM se arch for fields in the Lyman-c forest 
bv lOren & Wolfe (1995) found at most a marginal detection. 
If they exist, primordial magnetic fields may be detectable by 
Faraday rotation of the cosmic microwave background (CMB) 
polarisation durin g and shortly a fter the epoch of recombina- 
tion, as lKosowskv & Loebl 1119 96) proposed. 

IOhnoeta"fl J2002) proposed to use the CMB polarisation 
even for RM studies of nearby galaxy clusters. Should this 
speculative proposal become technically feasible, a lot of de- 
tailed information on intra-cluster magnetic fields could be ob- 
tained. 



1.3. Philosophy of the paper 

If magnetic fields are sampled in a sufficiently large volume, 
they can hopefully be regarded to be statistically homogeneous 
and statistically isotropic. This means that any statistical av- 
erage of a quantity depending on the magnetic fields does not 
depend on the exact location, shape, orientation and size of the 
used sampling volume. 

The quantity we are interested in this paper is the auto- 
correlation (or two-point-correlation) function (more exactly: 
tensor) of the magnetic fields. The information contained in 
the autocorrelation function is equivalent to the information 
stored in the power-spectrum, as stated by the Wiener- Khinchin 
Theorem (WKT). We therefore present two equivalent ap- 
proaches, one based in real space, and one based in Fourier 
space. The advantage of this redundancy is that some quanti- 
ties are easier accessible in one, and others in the other space. 
Further, this allows to crosscheck computer algorithms based 
on this work by comparing results gained by the different ap- 
proaches. 

The observable we can use to access the magnetic fields is 
Faraday rotation maps of extended polarised radio sources lo- 
cated behind a Faraday screen. Since an RM map shows basi- 
cally the line-of-sight projected magnetic field distribution, the 
RM autocorrelation function is mainly given by the projected 
magnetic field autocorrelation function. Therefore measuring 
the RM autocorrelation allows to measure the magnetic auto- 



1 The conventions to describe the spectra may differ in these articles 
from the one used here. 
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correlation, and thus provides a tool to estimate magnetic field 
strength and correlation length. 

The situation is a bit more complicated than described 
above, due to the vector nature of the magnetic fields. This im- 
plies that there is an autocorrelation tensor instead of a func- 
tion, which contains nine numbers corresponding to the cor- 
relations of the different magnetic components against each 
other, which in general can all be different. The RM autocor- 
relation function contains only information about one of these 
values, the autocorrelations of the magnetic field component 
parallel to the line-of-sight. However, in many instances the 
important symmetric part of the tensor can be reconstructed 
and using this information the magnetic field strength and cor- 
relation length can be obtained. This is possible due to three 
observations: 

1 . Magnetic isotropy: If the sampling volume is sufficiently 
large, so that the local anisotropic nature of magnetic field 
distributions is averaged out, the (volume averaged) mag- 
netic autocorrelation tensor is isotropic. This means, that 
the diagonal elements of the tensor are all the same, and 
that the off-diagonal elements are described by two num- 
bers, one giving their symmetric, and one giving their anti- 
symmetric (helical) contribution. 

2. Divergence-freeness of magnetic fields: The condition 
V-B = (B is the magnetic field) couples the diago- 
nal and off-diagonal components of the symmetric part of 
the autocorrelation tensor. Knowledge of one diagonal el- 
ement (e.g. from a RM measurement) therefore specifies 
fully the symmetric part of the tensor. The trace of the 
autocorrelation tensor, which can be called scalar mag- 
netic autocorrelation function w(r), contains all the infor- 
mation required to measure the average magnetic energy 
density Eg = w(0)/(8 ji) or the magnetic correlation length 
4 = £*w(r)M0). 

3. Unimportance of helicity: Although helicity is a crucial 
quantity for the dynamics of magnetic fields, it does not 
enter any estimate of the average magnetic energy density, 
or magnetic correlation length, because helicity only affects 
off-diagonal terms of the autocorrelation tensor. The named 
quantities depend only on the trace of the tensor and are 
therefore unaffected by helicity. One cannot measure helic- 
ity from a Faraday rotation map alone, since it requires the 
comparison of two different components of the magnetic 
fields, whereas the RM map contains information on only 
one component. 

In a realistic situation, the sampling volume is determined 
by the shape of the polarised radio emitter and the geometry 
of the Faraday screen, as given by the electron density and 
the magnetic field energy density profile. The sampling vol- 
ume can be described by a window function, through which an 
underlying virtually statistical homogeneous magnetic field is 
observed. The window function is zero outside the probed vol- 
ume, e.g. for locations which are not located in front of the ra- 
dio source. Inside the volume the window function scales with 
the electron density (known from X-ray observations), with the 
average magnetic energy profile (guessed from reasonable scal- 
ing relations, but testable within the approach), and - if wanted 



- with a noise reducing data weighting scheme. The effect of 
a finite window function is to smear out the power in Fourier- 
space. Since this is an unwanted effect one either has to find 
systems which provide a sufficiently large window or one has 
to account for this bias. Since the effect of a too small window 
on the results depends strongly on the shape and size of the 
window at hand, a detailed discussion of all the cases which 
can happen in practice is beyond the scope of this paper. It has 
to be done for each application at hand separately. Only ide- 
alised cases are discussed here for illustration. But generally 
one can state, that the analysis is sensitive to magnetic power 
on scales below a typical window size, and insensitive to scales 
above. 

The same magnetic power spectrum can have very different 
realizations, since all the phase information is lost in measuring 
the power spectrum (for a n instructive visualisation of this see 
Maron & Goldreich 2001). Since the presented approach relies 
on the power spectrum only, it is not important if the magnetic 
fields are highly organised in structures like flux-ropes, or mag- 
netic sheets, or if they are relatively featureless random-phase 
fields, as long as their power spectrum is the same. 

The autocorrelation analysis is fully applicable in all such 
situations, as long as the fields are sampled with sufficient 
statistics. The fact that this analysis is insensitive to differ- 
ent realizations of the same power spectrum indicates that the 
method is not able to extract all the information which may 
be in the map. Additional information is stored in higher order 
correlation functions, and such can in principle be used to make 
statements about whether the fields are ordered or purely ran- 
dom (chaotic). The information on the magnetic field strength 
(B 2 , which is the value at origin of the autocorrelation func- 
tion), and correlation length (an integral over the autocorrela- 
tion function) does only depend on the autocorrelation function 
and not on the higher order correlations. 

The presented analysis relies on having a statistically 
isotropic sample of magnetic fields, whereas MHD turbulence 
seems to be locally inhomogeneous, which means that small 
scale fluctuations are anisotropic with respect to the local mean 
field. However, whenever the observing window is much larger 
than the correlation length of the local mean field the auto- 
correlation tensor should be isotropic due to averaging over 
an isotropic distribution of locally anisotropic subvolumes. 
This works if not a preferred direction is superposed by other 
physics, e.g. a galaxy cluster wide orientation of field lines 
along a preferred axis. However, even this case can in prin- 
ciple be treated by co-adding the RM signal from a sample 
of clusters, for which a random distribution of such hypothet- 
ical axes can be assumed. In any case, it is likely that mag- 
netic anisotropy also manifests itself in the Faraday rotation 
maps, since the projection connecting magnetic field configu- 
rations and RM maps will conserve anisotropy in most cases, 
except alignments by chance of the direction of anisotropy and 
the line-of-sight. The presence of anisotropy can therefore be 
tested, which is discussed later in great detail. 

Since there are cases where already an inspection by eye 
seems to reveal the existence of magnetic structures like flux 
ropes or magnetic sheets, we briefly discuss their appearance in 
the autocorrelation and the area filling statistics of RM maps. 
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As already stated, the presence of such structures does not 
limit our analysis, as long as they are sufficiently sampled. 
Otherwise, one has to replace e.g. the isotropy assumption by a 
suitable generalisation. In many cases this will allow an analy- 
sis similar to the one proposed in this paper. We leave this for 
future work and applications where this might be required. A 
criteria to detect anisotropy statistically is given in this work. 

1.4. Structure of the paper 

In Sect. |2 the autocorrelation functions of magnetic fields and 
their RM maps are introduced, and their interrelation inves- 
tigated. In Sect. |3 the same is done in Fourier space, which 
has not only technical advantages, but also provides insight 
into phenomena such as turbulence. Faraday map signatures 
of magnetic structures like flux-ropes are briefly discussed in 
Sect. 0] Possible pitfalls due to observational artefacts are in- 
vestigated in Sect.|5] The conclusions in Sect.[6]summarise our 
main findings, and give references to the important results and 
formulae in detail. 

2. Real space formulation 

2.1. Basics 

The Faraday rotation for a line-of-sight parallel to the z-axis 
and displaced by x± from it, which starts at a polarised radio 
source at z s (jCj_) and ends at the observer located at infinity is 
given by 



RMOj.) = ci 



dz n t 



(x) B z (x) . 



(1) 



where ao = e 3 /(2nm 2 c 4 ), x = (x±,z), n e the electron density, 
and B the magnetic field strength. We assume in the following 
that any Faraday rotation due to a foreground as the Galaxy 
or the Earth's ionosphere is subtracted from the RM values, 
and only the RM of the Faraday screen is remaining. We also 
neglect any redshift effects, which can be included by inserting 
the factor (1 + Zredshift(z)) 2 into the integrand. 

The focus of this work is on the statistical expectation of 
the two-point, or autocorrelation function of Faraday rotation 
maps. This is defined by 



C RM (x^r ± ) = <RM(* ± )RM(* ± + rj>, 



(2) 



where the brackets indicate the expectation value of a statisti- 
cal ensemble average, which in practice may be replaced by a 
suitable average, e.g. over an observed RM map. For a given 
polarised background radio source of projected area Q. we de- 
fine the observable correlation function as 



— fdx 2 . 



RM(jCj.) RM(x ± + r ± ) , 



(3) 



where it is assumed that RM(xi) = for x± £ Q. An is taken 
here to be the radio source area, although other normalisations 
are imaginable, e.g. one might choose the area over which x±_ 
and x± + r ± are simultaneously within £1 This latter weight- 
ing would give a statistically unbiased estimator of Crm(/"j_), 
but we do not recommend its usage for the following reasons. 



The resulting estimator has a strong sensitivity to poorly sam- 
pled variances on scales comparable to the radio source diam- 
eter. On larger scales it is also ill-defined in the mathematical 
sense. The normalisation proposed here leads to results which 
do not suffer from this. It automatically down-weights the sig- 
nal from statistically insufficiently sampled baselines r±, and 
this bias can be controlled ("Sect. l3~5l . It further turns out that 
the Fourier-space formulation of RM statistics introduced in 
Sect.[3]requires the weighting scheme proposed here. 

The RM signal from different subvolumes of the Faraday 
screen will differ due to electron density and typical mag- 
netic field strength variations within the source. Such global 
variations can be regarded as variations of a window function 
fix), which mediates the relation between the observed RM 
signal and an underlying (rescaled) magnetic field, which is 
virtually homogeneous in a statistical sense. To be more spe- 
cific, we choose a typical position x K f within the screen (e.g. 
the centre of a galaxy cluster), and define n e o = n e (x K f) and 
Bq = (B 2 (x K f)) 1 ^ 2 . We then define the window function by 



/(*) = If* 



(4) 



where ^-{condition) is defined to be 1 if condition is true, otherwise 
0. g(x) = {B 2 (x)) l l 2 1 Bq is the dimensionless average magnetic 
field profile. In galaxy clusters a reasonable working assump- 
tion 2 for this ma ybejg(x) = (n < .(x) /n e f)) aB , for which we expect 
as ~ 1 (e.g. see Dolag et aLCOOll) . The function h(x ± ) allows 
us to assign different pixels in the map different weights 3 , e.g. 
in cases where the noise is a function of the position one might 
want to down-weight noisy regions 4 . If no weighting applies 
h(x±) = 1 everywhere. 

The expectation of the observed RM correlations are 



<C&(rJ> = ^ J>* jdr z 

f(x) fix + r) (B z (x)B z (x + r)), 



(5) 



with r = (r ± , r z ), and B(x) = B(x)/gix) the rescaled magnetic 
field. If properly rescaled, the average strength of the field is 
independent of the position. In that case, the rescaled magnetic 



2 We note that the average magnetic field strength profile can be 
tested a posteriori for consistency with the data, by comparing the 
observed values of RM 2 (jc) to the model expectations (see Sect. l2.4> . 

3 If hix ± ) j= 1 for some x ± e O the corresponding weight has to 
be introduced into the analysis. The most efficient way to do this is 
to make the data weighting virtually part of the measurement process, 
by writing RM(x ± ) — > RM w (x ± ) = RM(jc ± )/i(jc ± ) . For convenience, 
we drop the superscript™ again in the following, and just note that 
the analysis described here has to be applied to the weighted data 
RM w (x ± ). 

4 If <t(x ± ) is the noise of the RM map at position x ± a reasonable 
choice of a weighting function would be hix ± ) = cr /crix ± ). If the 
noise map itself could have errors, the danger of over-weighting noisy 
pixels with underestimated noise can be avoided by thresholding: If 
0" o is a threshold below which the noise is regarded to be tolerable, we 
propose to use hix±) = 1/(1 + cr^jj/cro). Another choice would be 
/j(x ± ) = l{o-(* x )«ro] which just cuts out regions which are recognised 
as too noisy. 
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field autocorrelation tensor should also be independent of posi- 
tion: 



field distribution, in the following called briefly isotropic tur- 
bulence): 



M i j{r) = {B i {x)Bj{x + r)). 



(6) 



w(r) = w(r) = 2M N (r) + M L (r) 



If the spatial variation of the window function is on much larger 
scales than the correlation length As of the magnetic fields, then 
Eq.|5]can be approximated to be 



1 d 

r 2 dr 



(r 3 M L (r)). 



(12) 



CrmOx) = <CRM( r -L)> = fl l C i( r i)' With fl l = a W e0 L > 



cm= r 

J -c 



dr z M zz {r), and r = (r ± , r z ). 



(7) 



In the last step Eq. [5] was used. Since the average magnetic 
energy density is given by (sb) = w(0)/(87r) the magnetic 
field strength can be determined by measuring the zero-point 
of w(r). This can be done by Faraday rotation measurements: 
The RM autocorrelation can be written as 



Here, we introduced the characteristic depth of the Faraday 
screen L - Vr/j/An, where V[f] - Jdx 3 f 2 (x) is the probed 
effective volume. We also introduced for convenience the nor- 
malised RM autocorrelation function C± which differs from 
Crm only by a geometry dependent factor a i . 

In the following we ignore the influence of the window 
function in the discussion, since for sufficiently large windows 
it only affects a \ . We therefore write B for B and keep in mind 
that our measured field strength Bo is estimated for a volume 
close to the reference location jc ref . At other locations, the av- 
erage magnetic energy density is given by Sb(x) = g 2 (x) B^. 
This approach assumes implicitly that typical length scales are 
the same throughout the Faraday screen. For sufficiently ex- 
tended screens, this assumption can be tested by comparing re- 
sults from different and separately analysed regions of the RM 
map. 

2.2. Isotropic magnetic correlation tensor 

The magnetic autocorrelation tensor for homogeneous 
isotropic turbulence, as assumed throughout the rest of this 
paper, can be written as 

My(r)=M N (r)£y + (M L (r) - M N (r)) ^ + M H (r) e ijk r k (8) 



CJrJ = ~ jdr t W( + r 2 ) = jc 



dr 



rw(r) 



(13) 



and is therefore just a line-of-sight projection of the magnetic 
autocorrelations. Thus the magnetic autocorrelations w(r) can 
be derived from C±(r±) by inverting an Abel integral equation: 



w(r) 



2 d r 
n r dr J r 

2 f 

Jr 



yC ± (y) 



C'Jy) 
yjy 2 — r 



(14) 



(15) 



where the prime denotes a derivative. For the second equation 
it was used that w(r) stays bounded for r — > oo. 

Now, an observational program to measure magnetic fields 
is obvious: From a high quality Faraday rotation map of a ho- 
mogeneous, (hopefully) isotropic medium of known geometry 
and electron density (e.g. derived from X-ray maps) the RM 
autocorrelation has to be calculated (Eq.0. From this an Abel 
integration (Eq.ll4lorll5> leads to the magnetic autocorrelation 
function, which gives (B 2 ) at its origin. Formally, 



(B 2 ) = w(0) 



2 r CJy) 
n Jo y 



(16) 



but this formulation is notoriously sensitive to noise. More sta- 
ble methods are presented later. 



(e.g.[Subramanian 1999) where the longitudinal, normal, and 
helical autocorrelation functions, Ml(t), M^(r), and Mu(r) re- 
spectively, only depend on the distance, not on the direction. 

The condition V B = leads to d I dr t Mijif) = (here and 2.3. Correlation volume, area, and length 
below we make use of the sum convention). This allows us to 
connect the non-helical correlation functions by 



M N (r)= l-^-(r 2 M L (r)) 
2r dr 



(9) 



For isotropic magnetic turbulence statements about integrals of 
w(r) and C±(r ± ) can be made. If correlations are short ranged, 
in the sense that M\Xr) r 3 — » for r — » oo, e.g. because of a 
finite size of the magnetised volume, then Eq.ll llimplies 



JSubramanianl 1 9991) . The zz-component of the magnetic auto- 
correlation tensor depends only on the longitudinal and normal 
correlations, and not on the helical part: 



V R = 



4n 



\drr 2 w(r) = Q. 
w(0) Jo 



(17) 



M„(r) = M h (r) 4 + Af N (r) -j with r - (r . r z ) , 



(10) 



which implies that Faraday rotation is insensitive to magnetic 
helicity. It is also useful to introduce the magnetic autocorrela- 
tion function 



This means that the magnetic autocorrelation volume Vb is zero 
and w(r) must therefore have positive and negative values if the 
magnetic field is non-zero. With the help of Eas.ll3lorll5lthis 
can be translated into a statement on the more directly measur- 
able RM autocorrelation: 



Mr) = (B(x) B(x + r)> = Miiir) , 



(11) 



IRM 



2n 
Cx(0) 



\dr ± 
Jo 



Cx(r x ) = 0. 



(18) 



which is the trace of the autocorrelation tensor, and depends 
only on r (in the case of a statistically isotropic magnetic 



Also C±(r) must have positive and negative values for a non- 
zero magnetic field, since the autocorrelation area A RM of the 
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RM map is zero. We note, that Eq.^]can conveniently be ex- 
pressed as 



^4rm — 



( fdx ± RM(jc ± )) 2 
Jdx ± RM^O 2 



(19) 



a form which should make clear that any foreground RM has to 
be removed before this quantity can satisfy Eq.^] Eqs. ll7llT%l 
and are a direct consequence of V-B = 0. These equations 
can be used to test how much a given observation (or model) 
deviates from giving a proper average of an isotropic ensemble. 

For example, the frequently used magnetic cell-model, in 
which cells of length-scale / ce n are filled with a from cell-to-cell 
randomly oriented but internally homogeneous magnetic field, 
does not fulfil V-B = 0. As a consequence it gives positive 
autocorrelation volumes and surfaces of the order Vb ~ ^ e]] 
and Arm ~ P ell respectively. Therefore it should be possible to 
exclude such an oversimplified model observationally. 

Other integral quantities of the turbulent magnetic field to 
look at are the existing non-zero magnetic and RM autocor- 
relation lengths. The magnetic autocorrelation length can be 
defined as 



A B 



= \dr 

%J -co 



w(r) C ± (0) 



— oo W(0) 



= 2 



w(Q) 



(20) 



where for the derivation of the last expression Eq. ^] or ^] 
can be used. Even in globally homogeneous turbulence, there is 
always a preferred direction defined by the local magnetic field. 
One can ask for the correlation length along and perpendicular 
to this locally defined direction and gets A\\ = jA B , and A ± = 
\A B so that A B = \{A\\ + 2A ± ). 5 

An observationally easily accessible length scale is the 
Faraday rotation autocorrelation length: 



tRM 



I dr ± 

%J — CO 



Xco 

Cx(0) K r drw(r) ' 



(21) 



From comparing Eqs.|20|and|^it is obvious that the RM cor- 
relation length-scale is not identical to the magnetic field auto- 
correlation length. As shown later, the RM correlation length 



5 The ratio of A\\/A ± = 2 seems to be small given the highly 
anisotropic nature of MHD turbulenc e (some recent progress on this 
can be found inlSridhar & Goldreicrll99ilGoldreich & SridhaJl997t 
iMaron & Goldreichl200lMCho et all2002l) . But it has to be noted, that 
there usually the magnetic fluctuations on top of a mean field are in- 
vestigated, whereas here the correlation lengths of an isotropically ori- 
ented magnetic field distribution sampled on length-scales above the 
scale of the locally present mean field of these studies, is investigated. 
The ratio A\\/Ax = 2 seems also to be in contradiction to the typical 
MHD turbulence picture of long and thin magnetic flux ropes with 
large aspect ratios. However, this is not necessarily the case. For a 
given position x ropc on a flux rope, the correlation length along and 
perpendicular will typically have a much larger ratio than Aj\/A ± = 2. 
But the position x mpc is a special position, since it was selected to be 
a high field strength region, whereas the correlation length Ay and A ± 
are defined by the statistical average over all positions in the volume. 
In order to make statements about the presence or absence of mag- 
netic filaments higher order statistics, beyond the two-point level used 
mostly in this work, have to be applied. 



is more strongly weighted towards the largest length-scales in 
the magnetic fluctuation spectrum than the magnetic correla- 
tion length. This is crucial, since in some cases these scales 
have been assumed to be identical, which could have led to 
systematic underestimates of magnetic field strengths. 

Having now defined two characteristic length-scales of the 
fields and their RM maps, suitable criteria testing the observed 
autocorrelation volume and area for statistical completeness 
can be formulated: |V° bs | « A\ and |AJ*£| « A 2 ^. We note 
that for strictly positive and therefore unphysical autocorrela- 
tion functions one would expect |V° bs | ~ A\ and |Aj*^| ~ A 2 ^ 
statistically. 

This can be turned around. For a sufficiently large RM map 
of a Faraday rotation screen with isotropic magnetic fields any 
homogeneous RM contribution from additional magnetised 
foregrounds can relatively accurately be measured. If there is a 
weak screen-intrinsic homogeneous magnetic field component, 
its z-component is 



Jd 2 x ± RM(x x ) 

«e0 Jd 3 xf(x ± ) 



(22) 



where in f(x) the average magnetic field profile may be set to 
g(x) = 1 everywhere for a truly homogeneous component. 

2.4. Testing the window function 

Now, all the necessary tools are introduced to test if the window 
function f(x) was based on a sensible model for the average 
magnetic energy density profile g 2 (x) and the proper geometry 
of the radio source within the Faraday screen z s {x±) (see Eq.|4}. 
Models can eventually be excluded a-posteriori on the basis of 



X 2 (*±) = 



RM(*J 2 



<RM(*x) 2 > ' 
where for the expected RM dispersion 



<RM(*J 2 > 



1 C 

-a 2 n 2 B 2 Q A B Jc 



dzf 2 (x) 



(23) 



(24) 



has to be used. As shown before Bo (e-g- by Eq. I16>. and As 
(e.g. by Eq. I20i can be derived for a given window function 



f(x) using CS^(rx). For a good choice of the window function, 



one gets 



1 



dx 2 L x 2 (x ± ) « 1 



(25) 



and larger values if the true and assumed models differ signif- 
icantly. The model discriminating power lies also in the spa- 
tial distribution ^ 2 (x ± ), and not only in its global average. If 
some large scale trends are apparent, e.g. that^ 2 (jcx) is system- 
atically higher in more central or more peripheral regions of 
the Faraday screen, then such a model for f(x) should be dis- 
favoured. This can be tested by averaging x 2 ( x ±) e -g- m radial 
bins for a roughly spherical screen, as a relaxed galaxy cluster 
should be, and checking for apparent trends. 

We note, that this method of model testin g can be re- 
garded as a refined Laing-Garrington effect dLaing l [1] 
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iGarrington et all 19 88): The more distant radio cocoon of a ra- 
dio galaxy in a galaxy cluster is usually more depolarised than 
the nearer radio cocoon due to the statistically larger Faraday 
depth. This is observed whenever the observational resolution 
is not able to resolve the RM structures. Here, we assume 
that the observational resolution is sufficient to resolve the RM 
structures, so that a different depth of some part of the radio 
source observed, or a different average magnetic energy profile 
leads to a different statistical Faraday depth (RM(jc_l) 2 ). Since 
this can be tested by suitable statistics, e.g. the simple^ 2 statis- 
tic proposed here, incorrect models can be identified. 

It may be hard in an individual case to disentangle the ef- 
fect of changing the total depth z s of the used polarised radio 
source if it is embedded in the Faraday screen, and the effect of 
changing B ( ° bs , since these two parameters can be quite degen- 
erate. However, there may be situations in which the geometry 
is sufficiently constrained because of additional knowledge of 
the source position, or statistical arguments can be used if a 
sufficiently large sample of similar systems were observed. 



normalised RM map is also identical to this, as a transformation 
of Eq.0shows, we can state 



(32) 



This Fourier-space version of Eq. ^] says, that the 2-d trans- 
formed RM map reveals the k z = plane of M zz (k), which in 
the isotropic case is all what is required to reconstruct the full 
magnetic autocorrelation w(k) = 2 C±(k). 

3.2. Power spectra 

We recall that the power spectrum P[F](k) of a function F(x) 
is given by the absolute-square of its Fourier transformation 
P[ F ](k) = \F(k)\ 2 . The WKT states that the Fourier transfor- 
mation of an autocorrelation function C[p](r), estimated within 
a window with volume V„ (as in Eq. [3}, gives the (windowed) 
power spectrum of this function, and vice versa: 



P [P] (k) = V n C [F] (k) . 



(33) 



3. Fourier space formulation 

3.1. Basics 

We use the following convention for the Fourier transformation 
of a «-dimensional function F(x): 



F(k) = Jd"xF(x)e ik - x 

F( X ) = -V fd n kP(k)e- ikx . 
(2n)» J 



(26) 
(27) 



The Fourier transformed isotropic magnetic autocorrelation 
tensor reads 

M u {k) = M N (k) |$ y - ^ J - is ijm ^ H(k) , (28) 

where we have directly used the V-B = condition in the form 
kj Mij(k) = to reduce the degrees of freedom to two compo- 
nents, a normal and helical part. The two corresponding spheri- 
cally symmetric functions in £-space are given in terms of their 
real space counterparts as: 

M N (k)= I d 3 rM N (r) e ik ' r =4n far r 2 M N (r) ^ r) (29) 



: Jd 3 rM N (r)e ik - r =47rjj. 

H(k)=^-M H (k)=^- \d\M H 
dk dkj 



kr 



(r)e ik - r , 



An C°° 

— \drr 2 M H {r) 
kJo 



kr cos(kr) - sin(fcr) 
kr ' 



(30) 



One can also introduce the Fourier transformed trace of the au- 
tocorrelation tensor w(k) = Mu(k) = 2 M^{k). A comparison 
with the transformed zz-component of the autocorrelation ten- 
sor 



M zz (k) = M N (k)(l -k 2 Jk 2 ) 



(31) 



reveals that in the k z = plane these two functions are identical 
(up to a constant factor 2). Since the 2-d Fourier transformed 



The WKT allows us to write the Fourier transformed auto- 
correlation tensor as 



1 



M u (k) = -(B i (k)B j (k)), 



(34) 



where V denotes the volume of the window function, which 
is for practical work with RM maps often the probed effective 
volume V = V[/j as defined in Sect. 12. II 

Thus, the 3-d magnetic power spectrum (the Fourier trans- 
formed magnetic autocorrelation function w{r)) can be directly 
connected to the one-dimensional magnetic energy spectrum in 
the case of isotropic turbulence: 



£#(£) dk — -^—r ~ - dk = n~rj dk . 



(2nY Sn 



2(2tt) 3 



(35) 



where we wrote w(k) = w(k) due to isotropy. The WKT also 
connects the 2-dimensional Fourier-transformed RM map with 
the Fourier transformed RM autocorrelation function: 



C x (Jtj.) = 



(\RM(kJ\ 2 ) 



(36) 



Thus, by comparing Eqs. [35] and [36] one finds that 
the magnetic energy spectrum is most easily measured from 
a given observation by simply Fourier transforming the map 
RM(*x), and averaging this over rings in £j_-space: 



«i A n (2. 



— f #|RM(*J| 2 
nf Jo 



(37) 



where k ± = k (cos ip, sin0). Eq. [Ogives a direct model in- 
dependent observational route to measure the turbulent energy 
spectrum. The average magnetic energy density can be easily 
obtained from this via 

pL ?W= f Ai y^, <38, 
Jo J a x A n (lnr 

where the last integration extends over the Fourier transformed 
RM map and can be done in practice by summing over pixels. 
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Fig. 1. Left: C±(r ± ) for power-law spectra of magnetic fluctuations w(k) = wo k~ s ~ 2 \\k\<k<k^) within the range k\ — 1 to k^ — 10 4 
with normalisation wq — 1 and different slopes s as labelled. For comparison, the thin line shows the RM correlation function 
expected for a single-scale power spectrum, as it can arise for isotropic linear force-free fields (C±(r ± ) - n (B 2 ) Jo(&f r ± )/(2 kp) 
with kp = 3, and (B 2 ) = 0.075; see Sect. l3.4> . Right: 1 - C±(r ± )/C±(0) in logarithmic units for the same power-law spectra. The 
thin lines are the asymptotic spectra given by the first two terms in Eg. 1441 



Also the correlation lengths can be expressed in terms of 

my. 



£°dkkw(k) 
§™ dkk 2 w(k) 



£dkwm 

" poo 

J dk k w(k) 



^RM = 2 J " (39) 

Jo """" Jo 
Thus, the RM correlation length has a much larger weight on 
the large-scale fluctuations than the magnetic correlation length 
has. Equating these two length-scales, as sometimes done in the 
literature, is at least questionable in the likely case of a broader 
turbulence spectrum. In typical situations (e.g. for a broad max- 
imum of the magnetic power spectrum as often found in hydro- 
dynamical turbulence) one expects Ab < Arm- Since the former 
is the one which enters the magnetic field estimates by using 
the measured RM-dispersion, 

2 <RM 2 > 



2C RM (0) 
fli Ab 



. , -,withL=-^, 

a^n L t LA B A n 



(40) 



using the easily measurable Arm instead of Ab likely underesti- 
mates the magnetic field strength. 

The isotropic magnetic autocorrelation function can be ex- 
pressed as 



w(r) = 



An 

(2^ 



o kr 



Similarly, the RM autocorrelation function can be written as 



CJr ± ) = -L 

47T 



dk k w(k) Jo(k r±) , 



(42) 



where J n (x) is the n-th Bessel function. In order to analyse the 
behaviour of the RM autocorrelations close to the origin, it is 
useful to rewrite the last equation as: 

kw(k) 



CM 



Jo 4n Jo 



kw(k) 



4n 



k(kr ± )). 



(43) 



The first term gives Cj_(0), and the second describes how 
C±(r±) approaches zero for r± — > oo. 



3.3. Power-law power spectra 

In many cases the small-scale magnetic energy spectrum is a 
power-law, say £g(£) = so& ' s (e-g- s — 5/3 for Kolmogorov- 
like turbulence, as expected if the magnetic fields were shaped 
by a mostly hydrodynamical turbulence) or wsik) = wok~ 2 ~ s 
for k\ < k < k2 (with k\ «c £2). For the behaviour of C±(r±) on 
such scales Eq.^]can be written as 



C ± (rJ = C x (0) - G(s) wo r[ + R M (r ± ) , 
where 

v(k) 



C x (0) 



r~k_m 

JO 47T 

2 S+2 n sTi 2 ^-) 



for < s < 2, and 



(44) 

(45) 
(46) 



= -1 fd. 
4ttJ 



dkk(w(k) - wo k- s - 2 ) (1 - J (Jtr)). (47) 



On small scales (more specifically for k\ <K l/r ± «; k%) and 
for well behaved power spectra outside k\ <c k <K k-i the term 
R[w](r) is negligibly small. We therefore propose to fit 



(41) CRM(r±) ~ Co - C\ r s x 



(48) 



to the inner part of an observationally determined RM cor- 
relation function. From this, the turbulence spectral index 
s, Cj_(0) = Co/ai, and the power-law normalisation wq = 
C\l{a\ G(s)) can be inferred. A rough estimate of k\ and ki, 
the scales on which the spectrum deviates from the power-law, 
can also be obtained from finding the r ± -values, where the fit 
becomes poor. A more accurate determination of these scales 
can always be done in the Fourier domain (see Eas. 1321 l36l and 
I37> . We therefore recommend the usage of Eq. |48]more for 
consistency checks and rapid and rough fit-by-eye diagnostics 
of the steepness of the magnetic power spectrum rather than for 
high-precession analysis. 
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The shape of the RM correlation function close to the ori- 
gin allows a direct read-off of the type of magnetic turbulence. 
A top-down scenario, where most of the energy resides on 
large scales (s > 1) and the smaller scales are populated by 
a turbulent cascade as in the Kolmogorov-, Kraichnan-, and 
Goldreich-Sridhar-phenomenologies, leads to a fiat cusp at the 
origin, and a convex shape near to it. A bottom-up magnetic 
turbulence scenario, where the fields originate on small scales 
and are enlarged by shear flows or other inverse cascade ac- 
tions (s < 1), leads to a sharp cusp at the origin, and a concave 
slope next to it. A spectral energy distribution with as much en- 
ergy on small as on large scales (s = 1) leads to a linear cusp 
at the origin. The behaviour of C±{r ± ) for these three cases is 
illustrated in Fig.^ 

3.4. Helical correlations & force-free fields 

Faraday rotation maps do not contain information about the he- 
lical part of the autocorrelation tensor. Therefore, additional in- 
formation is required in order to be able to measure the helical 
correlations. For example any relation between the helical and 
non-helical components would be sufficient. 

In order to give an example for such additional informa- 
tion, we discuss the case of force-free fields (FFFs). The con- 
dition for FFFs reads VxB = k F B , where k F can in general 
be a function of position. For simplicity, we restrict k F to be 
spatially constant. Such so-called linear FFFs lead to a very 
simple structure of the components of Mjj(k): For k + k F all 
components vanish, and for k = k F one gets M^ike) - -H(k F ), 
leaving the magnetic energy density (or the helicity) as the only 
remaining free parameter for a given characteristic wave-vector 
k F . FFFs are therefore also called maximally helical fields. 

From the fact that for a linear FFF only one spherical shell 
in wave-vector space is populated with magnetic power, the 
spatial autocorrelation function is easily obtained as 

w F (r) = (B 2 )sm(k F r)/(k F r) (49) 

(Subramanian 1999) and the RM autocorrelation as 

CJrJ = n (B 2 ) ] (k F r0/(2 k F ). (50) 

This function is shown on the left side of Fig.^ 

3.5. Finite window functions 

Here, we discuss the effect of a finite window function on mag- 
netic field estimates, in order to possibly correct for the bias 
made with the robust weighting scheme introduced in Sect. 12. II 
Taking a finite window function f(r) into account, Eq. I32lbe- 
comes 



Cf(*J = i jd 3 qw{q)^ W{k ± 
where we introduced 



q) 



(51) 



(52) 



and used the identity M zz {k) - \w{k) (1 - k 2 /k 2 ). Without the 
term q 2 ± /q 2 the convolution integral in Eq.|^ would describe 
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Fig. 2. Response in C±(k±) to a single-scale magnetic power 
signal w(k) = 5{k - p) at wavenumber p for different radii R of 
a cylindrical window function. 



a redistribution of the magnetic power within Fourier space, 
which conserves the total magnetic energy. But the term q^/q 2 
leads to some loss of magnetic power. 

In any situation in which there is substantial magnetic 
power on scales comparable or larger than the window size 
Eq. |^can be used to estimate the response of the observa- 
tion to the magnetic power on a given scale p by inserting 
w(q) = 6{q-p). Ideally, this is then used within a matched-filter 
analysis or as the response matrix in a maximum-likelihood re- 
construction of the underlying power-spectra. The computation 
of the response matrix relating input power w{q) and measured 
signal C° bs (£j_) can be cumbersome since in general a 2- or 
3-dimensional integral (depending if one uses the delta func- 
tion) has to be evaluated for each matrix element. Therefore 
we restrict our discussion here to three highly symmetric, ide- 
alised cases, and an approximative treatment of a more realistic 
configuration, which should give a feeling for the general be- 
haviour. 

A cylindrical window: Suppose a circular radio source 
with radius R is seen through a very deep Faraday screen, so 
that the depth L z can be approximated to be infinite long. The 
window function f(x) = l{ Xx <R] leads to 



W(k) = 



J?(*x*) 

Tlk\ 



S(k z ), 



(53) 



where 5(k) is the Dirac's delta function. Inserting this into 
Eq.l5T1 gives 



P r 2n J K >M + p 

2ttJ 



■2pk ± cos (pR) 



k]_ + p 2 - 2 p k ± cos < 



(54) 



which is shown in Fig.|2 One can clearly see that with increas- 
ing window size the response becomes more and more delta- 
function-like. 

A sheet-like window seen edge-on: Suppose a very elon- 
gated radio source is seen through a deep Faraday screen, so 
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Fig. 3. Response in C±(k± (cos <p, sin (/>)) to a single-scale mag- 
netic power signal w(k) = 6(k-p) at wavenumber p for a sheet- 
like window seen edge-on with diameter L x = 10 /p for differ- 
ent wavevector orientations <p with respect to the sheet normal. 

that the window function is approximated by f(x) = l(o<. t <L v |- 
This gives 

W(k) = A Lx (k x )8(ky)6(k z ), (55) 

where we introduced a shortcut for the Fourier transformed 1- 
d-box-window of size L: 

'"id / / 'i 

6(k) for L -> oo . (56) 



L sin z (fcL/2) 

Ai (k) = — 

In (kL/2) 2 
This leads to a response of 



l{ky< P } ( A L x (k X 



t) + A L Jk x + t)) 



(57) 



where we wrote t = (p 2 - k 2 ) 1 ^ 2 for brevity. The anisotropic 
response of this window is much more delta-function like for 
wavevectors oriented along the sheet axis than oriented perpen- 
dicular to it, as can be seen in Fig. [3] 

A sheet-like window seen face-on: Suppose the window 
is an infinitely extended homogeneous layer of thickness L z in 
z-direction, e.g. a magnetised skin layer of a large radio source, 
so that f(x) = 1{o< z <l ; |, and 

W(k) = A L: (k z )6 2 (kJ, (58) 

where 6 2 (kj_) is the 2-d Dirac's delta function. 
The response function 

k\ 



C°j_ s (k x ) - l{k ± <p] ■ 



■k 2 J 



(59) 



is isotropic, and is shown in Fig.@] In this geometry, power gets 
scattered out of the k z = plane due to the finite window size 
in z-direction. This leads to a loss of magnetic power in any 
measurement, which does not correct for this bias. 

The observed magnetic energy density estimated with the 
help of Eq.^JOcan be shown to be related to the real magnetic 
power spectrum via 



sf s = dks B (k)D(k), 



(60) 
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Fig. 4. Response in C±(k ± ) to a single-scale magnetic power 
signal w(k) = 6{k-p) at wavenumber p for a sheet-like window 
seen face-on for different diameters L 7 . 



where the weighting function 
3 



D(k) = 



2kL, 



iF 2 
kL 7 



1 1 

2' 2' 



2; - 







f) 


- 1 







V(16/3) 2 + (fcL z ) 2 



kL^ 

~Y 

3kL z /\6 ; kh, <K 1 
1 ; kL~ » 1 



(61) 



(62) 



describes the relative contribution of different parts of the spec- 
trum. Note that ssik) D(k) is not the observed power spectrum, 
but only the contribution to the magnetic energy estimate, that 
a Ft, denotes the hypergeometric function, and that the asymp- 
totic approximation has an overall accuracy of better than 5%. 
Since D(k) < 1 everywhere, it is obvious that the derived mag- 
netic fields are underestimated, especially if there is substantial 
magnetic power on scales comparable and larger than the win- 
dow size. 

In principle, it is possible to correct for any bias, if the win- 
dow function is reliably known and if the statistical sampling 
is sufficiently good even on the large scales so that dividing the 
observed magnetic power spectrum by the weighting function 
gives sensible results (and not just amplifies the noise). 

Approximative treatment of a realistic window: In a re- 
alistic situation, often a relatively small sized radio galaxy is 
seen through a deep Faraday screen. In such a case the depth 
can again be approximated to be infinite for the purpose of the 
Fourier-space window: 



W(k)=d(k z )W ± (ki.), 

where we introduced the projected Fourier-window 

lA(*x)i 2 



W ± (k±) = 



(2^) 2 A [/J 



which results from a projected window function: 
/ 2 (* ± ) = ^dzf(x)IL z Wx\hA [fA = jd 2 x ± fl(x ± ) . 



(63) 



(64) 



(65) 



T. A. EnBlin and C. Vogt: Faraday Rotation Screens 



11 



Here, L z is an arbitrary but fixed reference length, e.g. the typ- 
ical source size L z = V[/]/Aq. In the case that the observation- 
ally measured magnetic power spectrum results from a spheri- 
cal average over the data 



1 



2k 



w ob \k) = — \d(f>2 Cf s {k (cos 0, sin (/>)) 
2 7i Jo 



(66) 



the response to a delta function-like magnetic power-spectrum 
w(q) = 6(q - p) is given by 

Wj_(?j.) M\k ± -p\<q x <k±+p) 



wf\k ± ) 



2 p r,i 

— Id q ± 



2 -l fl 



*j- + p q jM W ± (q(co&(f), sin^)) 
dq -. 

"'I y]4q 2 p 2 -(q 2 +p z -k 2 J 2 



(67) 



(68) 



This approximative response can be easily computed numeri- 
cally for any model window function. In many cases it should 
be sufficiently accurate to estimate the effect of a finite ob- 
servational window on the derived magnetic power spectrum. 
As a consistency check, we verified that the limit of an in- 
finitely extended radio source, which can be written as u — » 
in W ± (q x ) = 1( 9± <„|/(7tm 2 ), leads to wf s (q) = 5{q - p) as it 
should. 



3.6. Testing isotropy 

Since isotropy of the magnetic field statistics is a crucial ingre- 
dient of the proposed analysis, it is important to test if indica- 
tions of anisotropy are present, and to see how anisotropy can 
affect the results. 

Anisotropy can manifest itself in two different ways: 
(a) The Fourier space magnetic power distribution can be 
anisotropic, by being not only a function of k but a full 
function of k, and (b) the magnetic power tensor itself can 
be anisotropic. Certainly both flavours of anisotropy can be 
present simultaneously. However, their effects can be well sep- 
arated, so that we discuss them one by one. 

Before doing so, we note that the relation 



Ci.(ifcx) = M zz (k ± ,0) 



(69) 



is completely independent on assumptions on isotropy. The 
RM power spectrum still reveals the k z = plane of the zz- 
component of the magnetic power tensor. The condition of 
isotropy had allowed to use the measured information as a rep- 
resentative probe of (a) the full Fourier-space and (b) the other 
diagonal elements of the magnetic tensor. 

(a) Anisotropic power spectrum: Eq. |69] shows that an 
anisotropic power spectrum can be detected, since it leads very 
likely to an anisotropic RM power map if the anisotropy is not 
aligned with the z-direction by chance. Since the latter can not 
be excluded, it is hard to prove isotropy. On the other hand, a 
perfect alignment of the line-of-sight and the anisotropy axis is 
not very likely. By studying a number of independent Faraday 
screens, an anisotropic power spectrum can be ruled out on a 
statistical basis. Furthermore, by co-adding the signals of sev- 
eral systems, statistical isotropy can be enforced, even if an in- 
dividual system is anisotropic. However, in order to be able to 



co-add different observations, the window functions have to be 
well understood. Especially the scaling of the average magnetic 
field energy density with location within the Faraday screen 
and from screen-to-screen should be known. Since this is still 
poorly known, it is worth to check if indications of anisotropy 
are present in every dataset itself. 

A good way to check for indications of anisotropy is by 
eye inspection of maps of C ± (k±) or equivalently C±(r ± ) or by 
comparing profiles which were calculated using different angu- 
lar slices. A more quantitative estimate of apparent anisotropy 
can be obtained by the use of multipole moments. Since the 
dipole moment vanishes due to mirror symmetries in C±{k ± ) 
and C±(r ± ), the first non-trivial multipole is the quadrupole 
moment: 



e: 



!?-/■ 



d k± C ± (k±) (2 kt kj - k\ tfy) k ± y 



(70) 



(and similarly for C±(r ± )). Here we allowed for a weighting 
factor k/ in order to balance the contributions from different 
scales. 2 <y) This should be compared to the second moment of 
the (weighted) distribution 



P (y) 

r 2 



d z k ± C ± (kJ k]_ 7 (and similarly for CJrJ), (71) 



e.g. by calculating the ratio R (y) = (\Q<g\ + \Q ( $\)/(2 P (7) ). 
In an isotropic case this number should be close to zero, in 
a strongly anisotropic case it can become comparable with 
one. We suggest to apply this test to the real-space and the 
Fourier-space data, since in the last case isotropy is mostly 
tested on large spatial scales, and in the second case on small 
spatial scale. Attention has to be given to the fact that the 
real-space quadrupole moment is sensitive (for small y) to any 
non-circularity of the window function which can affect large 
r x . The Fourier-space quadrupole moment can be affected (for 
small y) by any ellipticity of the synthesised beams of the ob- 
servations since this manifests itself at large k ± . Therefore the 
integration range for the quadrupole moments and the second 
moments may be better restricted to intermediate radii in real- 
and in Fourier-space. 

(b) Anisotropic tensor: In the anisotropic case, the only 
constraint on the magnetic autocorrelation tensor is k t Mifk) = 
due to the divergence-freeness of the magnetic fields. This 
translates for the z-components into k x M xz (k) + k y M yz (k) + 
k z M zz (k) = 0, which leaves the observable M zz (k ± ,0) abso- 
lutely unconstrained since k z = 0. M zz can therefore be an arbi- 
trary function of k ± . However, if it is not circularly symmetric, 
this can be detected with the methods described above. 

In order to have a working example of an anisotropic part 
Mfj of the magnetic tensor we assume that a preferred direction 
a(k) exists, so that 

ki ai k t 



Mf j (k) = \a i - 



k 2 



kj ai k\ 
~1 2 ~ 



(72) 



This is an intrinsically anisotropic tensor, which fulfils 
0. The Faraday-observable component is 
a 2 (k ± ). Its influence on the RM statistics can 
not be discriminated from an isotropic contribution if a{k) de- 
pends only on \k\. The assumption of isotropy would therefore 



kiMf.(k) 
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lead to an incorrect estimate of the field strength, since the 
measured zz-component is assumed to be representative for all 
components. However, if the signal from a number of similarly 
anisotropic Faraday screens are co-added the errors compen- 
sate statistically, if no correlation of the anisotropic direction 
and the line-of-sight are present. 

Furthermore, since any anisotropy of the magnetic power 
tensor should have a physical cause, e.g. a large-scale 
shear flow in the Faraday active medium, an accompanying 
anisotropic power spectrum is very likely, which can princi- 
pally be detected by the methods described above (a). 

Finally, if anisotropy turns out to be inherently present in 
Faraday screens, one might replace Eas. ISland I28lbv a more 
complex, anisotropic model in order to be able to extract infor- 
mation from individual screens. In that case this work may help 
as a guideline for such a more elaborate analysis. 

4. Magnetic structures 

4.1. Autocorrelation 

The possibility exists that the magnetic fields of a Faraday 
screen consist of several distinct magnetic structures like flux 
ropes, magnetic tori etc. If the positions and orientations of the 
structures can be regarded as statistically independent the mag- 
netic autocorrelation function can be written as 



(73) 



where all types s of structures present with space density n s are 
summed up. A structure s with field configuration B s (x) has an 
intrinsic isotropically averaged (unnormalised) autocorrelation 
function 



W s (r) 



d 3 xB s (x)-B s (x + flr), 



(74) 



where the first integration covers the unit sphere. 

For a magnetic structure, which consists of a mostly con- 
stant magnetic field B s within the volume V s , and negligible 
field strength elsewhere, the autocorrelation function is asymp- 
totically for small r 



W s (r) = B 2 s V s (l-r/l s ), 



(75) 



where l s is a typical length-scale of the structure, roughly given 
by l s ~ V s /A s with A s the surface area of the structure. If only 
a single type of structure is present, we get asymptotically 



w(r) = B 2 s Ti B (l-r/l s ), 



(76) 



where tjb = n s V s is the magnetic volume filling factor. 

In order to calculate the RM autocorrelation of such a 
Faraday screen, we use as a toy model w(r) from Eq. ^6] as 
long as r < r max , and otherwise w(r) = 0. Eq.^] would then 
requires r max = 4 l s , but the actual choice is only important for 
numerical values of constants of proportionality, and not for 
the qualitative shape of the RM autocorrelation function at the 
origin. Integrating Eq.^]leads to an asymptotic expansion of 
the form 

\2 



C±(r ± ) = C 



Cl+ C 2 ln(-) 
r± 



(77) 



which gives a flat central slope, nearly a parabola, but still hav- 
ing a tiny logarithmic cusp. The constants Co, C\, and C2 de- 
pend on the details of the outer slope of w(r), e.g. on the ratio 
r max/h> so m at their numerical values are model dependent. 

We summarise that a Faraday screen built from structural 
elements with internally constant magnetic fields, and only a 
single characteristic length-scale leads to a flat central autocor- 
relation function, with at most a logarithmic cusp of the form 
given by Eq.I77I 

4.2. Filling factor 

Although there exist characteristic shapes of the RM autocor- 
relation function C x (r±) in the case of a patchy magnetised 
Faraday screen, as demonstrated in Sect l4.ll the presence of 
such patches can not be deduced from C±(r±) alone. Since the 
phase information is missing, the special form of the cusp aris- 
ing from magnetic structures as given by Eq. El can not t> e 
distinguished from a complete random phase turbulence with 
steep power-law like spectra with spectral index s m 2, as can 
be seen from comparison with Eq.l48l 

In order to measure the patchin ess of the magnetic field 
distribution in galaxy clusters Cla rke et al.l (12001 ) used the area 
filling factor t]rm of the line-of-sight of extended radio sources 
which do not show any RM due to the Faraday screen. For their 
sources, they concluded that 77rm > 95%. 

If the magnetic fields are in flux-rope like structures, with 
typical length l\\ and diameter l ± , the cross section of a flux-rope 
to be intersected by a line-of-sight is of the order l\\ l x . Their 
volume filling factor is tjb ~ l\\ l]_ n mpe . If their locations can be 
regarded as being uncorrelated, the number K of flux ropes in- 
tersected by a line-of-sight of length Li os is Poisson-distributed: 
P(K) = A* exp(-A)/(A"!) with A * Tj B Li os /l h From that it fol- 
lows by inserting K = that 



A ln(— i 

Lios 1 - ?7rm 



(78) 



For filaments of length l\\ a 10 kpc and lines-of-sight of 
Lios ~ 500 kpc only a small subvolume t/b > 0.05 of the clusters 
actually needs to be magneti sed in order to give the large area 
filling factor found bv lClarke et ail fcOOlh . 

Another constraint for the magnetic filling factor can be 
obtained from energetic arguments. The magnetic field energy 
density in magnetised regions can be expected to be below the 
environmental thermal energy density e±, since otherwise a 
magnetic structure would expand until it reaches pressure equi- 
librium. Since the autocorrelation analysis of RM maps is able 
to provide the volume averaged magnetic field energy density 
(sb), the magnetic volume filling factor can be constrained to 
be 



(eg) 



0.8 ■ 10" 



(B 2 ) 
juG 2 



n e kT 
10- 3 cm- 3 keV 



(79) 



In cases of relatively strong average magnetic energy densities, 
as e.g. in cooling flow clusters, or in case of physical arguments 
requiring a much lower than equipartition field strength, this 
can give a tight constraint on the volume filling factor. 
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5. Observational artefacts 

5.1. Beam smearing 

The finite size Zb ea m of a synthesised beam of a radio interfer- 
ometer should smear out RM structures below the beam size, 
and therefore can lead to a smooth behaviour of the measured 
RM autocorrelation function at the origin, even if the true au- 
tocorrelation function has a cusp there. Substantial changes of 
the RM on the scale of the beam can lead to beam depolar- 
isation, due to the dif ferentially rotated polari s ation vectors 
with in the beam area dConwav & Stromlll985t lLaind ll988: 
iGarrington et alJlT988h . Since beam depolarisation is in prin- 
cipal detectable by its frequency dependence, the pre sence of 
sub-beam structure can be noticed , even if not resolved ( Tribble 
Il99ll iMelrose & Macauartlll998l) . The magnetic power spec- 
trum derived from a beam smeared RM map should cut-off at 
large k ~ 7T/l beam . 

5.2. Noise 

Instrumental noise can be correlated on several scales, since ra- 
dio interferometers sample the sky in Fourier space, where each 
antenna pair baseline measures a different fc ± -vector. It is diffi- 
cult to understand to which extend noise on a telescope antenna 
baseline pair will produce correlated noise in the RM map, 
since several independent polarisation maps at different fre- 
quencies are combined in the map making process. We there- 
fore discuss only the case of spatially uncorrected noise, as it 
may result from a pixel-by-pixel RM fitting routine. This adds 
to the RM autocorrelation function 

C^Jx ± , rj = ai^ixj 6 2 (r ± ) . (80) 
In Fourier space, this leads to a constant error for w(k) 

HWse(k) = 2 <Cr| M noise > (81) 

and therefore to an artificial component in the magnetic power 
spectrum s°g S (k) increasing by k 2 . 

If it turns out that for an RM map with an inhomogeneous 
noise map (if provided by an RM map construction software) 
the noise affects the small-scale power spectrum too severely, 
one can try to reduce this by down- weighting noisy regions 
with a suitable choice of the data weighting function h(x ± ) 
which was introduced in Sect. l2.1l for this purpose. 

5.3. RM steps due to the nn -ambiguity 

An RM map is often derived by fitting the wavelength-square 
behaviour of the measured polarisation angles. Since the po- 
larisation angle is only determined up to an ambiguity of nn 
(where n is an integer), there is the risk of getting a fitted RM 
value which is off by m ARM from the true one. m is an integer, 
and ARM = n(A 2 . - /l^, ax ) _1 is a constant depending on the 
used wavelength range from A m [ n to A max . 

This can lead to artifical jumps in RM maps, which will af- 
fect the RM autocorrelation function and therefore any derived 
magnetic power spectrum. In order to get a feeling for this we 



model the possible error by an additional component in the de- 
rived RM map: 

RM amb (jc ± ) = J] tm ARM l {x±eCli] , (82) 

i 

where Q,- is the area of the z'-th RM patch, and m, is an integer, 
mostly +1 or -1. Assuming that different patches are uncorre- 
lated, the measured RM autocorrelation function is changed by 
an additional component, which should be asymptotically for 
small r x 

C£*(rJ = ARM 2 77 amb f 1 - , (83) 

where ^ am b is the area filling factor of the ambiguity patches in 
the RM map, and / am b a typical patch diameter. 

Comparing this with Eq.|48] shows that the artificial power 
induced by the «7r-ambiguity mimics a turbulence energy spec- 
trum with slope 5 = 1, which would have equal power on all 
scales. A steep magnetic power spectrum can therefore possi- 
bly be masked by such artifacts. 

Fortunately, for a given observation the value of ARM is 
known and one can search an RM map for the occurrence of 
steps by ARM over a short distance (not necessarily one pixel) 
in order to detect such artifacts. 

6. Conclusions 

We have investigated the statistics of Faraday rotation maps on 
the level of the autocorrelation function and the power spec- 
trum. We proposed ways to analyse extended Faraday maps in 
order to reconstruct the magnetic autocorrelation tensor (Eqs.[6] 
and I34t from which quantities like the average field strength, 
the magnetic energy spectrum, and their autocorrelation length 
can be obtained (Eq. I20> . We showed that under the assump- 
tion of isotropy of the observed magnetic field ensemble the 
symmetric part of the magnetic autocorrelation tensor (Eqs.|8] 
and 12 8 \ can be reconstructed. This makes use of the condition 
VvB = and the additional assumption (which can be tested a- 
posteriori) that the gradient scale of the electron density (e.g. 
the core radius of a galaxy cluster) is much larger than the 
typical field length-scale. The anti-symmetric or helical part 
of the magnetic correlation tensor can only be measured if ad- 
ditional information is available, e.g. in the case of force-free 
fields (Sect.E3). 

The assumption of isotropy of the magnetic field statistics 
should be justified in cases where a sufficiently large volume 
of the screen is probed. In principle, it can also be tested by 
searching for non-circular distortions of the 2-dimensional au- 
tocorrelation function (Sect. l3~oY 

A further test for statistical isotropy and sufficient sampling 
of the field fluctuations is the fact that if these conditions are 
given in a finite Faraday screen (which cannot maintain in- 
finitely long correlations) the rotation measure (RM) autocor- 
relation area (Eq. 1181 1 has to vanish. This means that there is a 
balance between the positively and negatively valued areas of 
the autocorrelation function. In practice, one would require it to 
be much smaller than the RM autocorrelation length squared. 
We note that e.g. the popular magnetic cell-model, in which 
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cells are filled by from cell-to-cell independently oriented and 
internally homogeneous magnetic fields, does not have these 
properties, since it violates the required V\B = condition. 

Our approach is meant to be applied directly to real data. 
Effects of incomplete information, due to the limited extent of 
polarised radio sources, are properly treated in form of a win- 
dow function (Eqs.@J[5] and I51ll. This window function con- 
tains additional information on the screen geometry, and allows 
for proper bookkeeping of data weighting, in case of noisy data 
being analysed. Since the window function also requires some 
working hypothesis about the average magnetic energy density 
profile of the Faraday screen, ways to test it a posteriori are 
sketched (Sect. 12.41. 

The most efficient way to analyse Faraday rotation maps 
leads through Fourier space (Sect. [3}- The Fourier transforma- 
tion of a map gives direct insight into the magnetic energy spec- 
trum (Eas. 1351 and I37>. which fully specifies the magnetic au- 
tocorrelation function (Ea.l32>. Many quantities of interest can 
be obtained from it, such as the average field strength (Ea.l38>. 
the correlation lengths (Eq.l39L and the bias resulting from the 
used window function (Sect. l3~5Y 

The Fourier domain formulation gives also important in- 
sight into the real space behaviour of the RM autocorrelation 
function (Sect. 13 . 3I >: A power-law magnetic energy spectrum 
leads to a cusp at the origin of this function, where the shape 
of the cusp is determined by the power law index s. A steep 
spectrum s > 1 (as e.g. expected for turbulent cascades) leads 
to a flat cusp, whereas a flat spectrum s < 1 gives a pronounced 
sharp cusp. The limiting case s — 1 with equal power on all 
scales leads to a linear (decreasing) behaviour of the RM auto- 
correlation function close to the origin. 

Such cusps of the RM autocorrelation (or power-law spec- 
tra in Fourier space) can be signatures of turbulent cascades, but 
they can also occur in other situations. We demonstrated that a 
Faraday screen which is composed of finite magnetic structures 
of roughly constant field strength lead to a flat cusp, too (Sect. 
14. 1> . We show ways to constrain the magnetic volume filling 
factor for such magnetic field models (Sect. 14.2V Observational 
artifacts, like noise or jumps in the measured RM values due 
to the so called «7r-ambiguity, are able to produce sharp cusps 
(Sect.|5}- We therefore stress the importance to check maps for 
such distortions and describe ways to do this. 

A very important result of this work is that the magnetic 
autocorrelation length Ab (Eq.I20> is in general not identical to 
the autocorrelation length of the RM fluctuations Arm (Eq. l21> . 
/Irm is much more strongly weighted towards the large-scale 
part of the magnetic power spectrum than Ab (Eq.l39>. In typi- 
cal astrophysical situations a broad spectral energy distribution 
can be expected so that Ab will be much smaller than Arm- 
Since As enters the classical formulae to estimate magnetic 
field strength from Faraday measurements (B 2 ) oc (RM 2 ) /Ab 
(see Ea.l40>. but /Irm is sometimes used instead, we expect that 
several published magnetic field estimates from Faraday rota- 
tion maps are actually underestimates. 

We hope that our work aids and stimulates further obser- 
vational and theoretical work on the exciting field of Faraday 
rotation measurements of cosmic magnetic fields in order to 



give us deeper insight in their fascinating origins and roles in 
the Universe. 
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